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We study the process e + e _ — > tv + tt~ J/ip at a center-of-mass energy of 4.260 GeV using a 525 pb _1 data 
sample collected with the BESIII detector operating at the Beijing Electron Positron Collider. The Bom cross 
section is measured to be (62.9 ± 1.9 ± 3.7) pb, consistent with the production of the Y (4260). We observe a 
structure at around 3.9 GeV/c 2 in the n ± J/i/j mass spectrum, which we refer to as the Z c (3900). If interpreted 
as a new particle, it is unusual in that it carries an electric charge and couples to charmonium. A fit to the ir^ J/tp 
invariant mass spectrum, neglecting interference, results in a mass of (3899.0 ±3.6 ±4.9) MeV/c 2 and a width 
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of (46 ± 10 ± 20) MeV. Its production ratio is measured to be R = CT ( e + e - -»7r+7r- JMj 

(21.5 ± 3.3 ± 7.5)%. In all measurements the first errors are statistical and the second are systematic. 



PACS numbers: 14.40.Rt, 14.40.Pq, 13.66.Bc 



Since its discovery in the initial-state-radiation (ISR) pro- 
cess e + e~ — >■ 7isR7r + 7r - J/ip J]]], and despite its subsequent 
observations |2j-|5|], the nature of the y(4260) state has re- 
mained a mystery. Unlike other charmonium states with the 
same quantum numbers and in the same mass region, such 
as the V(4040), ^(4160), and ^(4415), the y(4260) state 
does not have a natural place within the quark model of 
charmonium [6]. Furthermore, while being well above the 
DD threshold, the F(4260) shows strong coupling to the 
7r + 7r~ J I tp final state 131, but relatively small coupling to open 
charm decay modes 18I 412I1 . These properties perhaps indicate 
that the F(4260) state is not a conventional state of charmo- 
nium lfl3ll . 

A similar situation has recently become apparent in the bot- 
tomonium system above the BB threshold, where there are in- 
dications of anomalously large couplings between the T(5S) 
state (or perhaps an unconventional bottomonium state with 
similar mass, the F b (10890)) and the 7r+7r-T(15, 2S, 35) 
and 7r + 7r"/i 6 (lP,2P) final states Q Q. More sur- 
prisingly, substructure in these ir + 7r~T(lS,2S,3S) and 
7r + 7r~ hb(lP, 2P) decays indicates the possible existence of 
charged bottomoniumlike states 11611 . which must have at least 
four constituent quarks to have a non-zero electric charge, 
rather than the two in a conventional meson. By analogy, 
this suggests there may exist interesting substructure in the 
y(4260) — > 7r + 7r~ J/tp process in the charmonium region. 

In this Letter, we present a study of the process e + e~ — ^ 
7r + 7r~ J/tpata center-of-mass (CM) energy of s/s = (4.260± 
0.001) GeV, which corresponds to the peak of the F(4260) 
cross section. We observe a charged structure in the 7r ± J / tp 
invariant mass spectrum, which we refer to as the Z c (3900). 
The analysis is performed with a 525 pb _1 data sample col- 
lected with the^BESIII detector, which is described in de- 
tail in Ref. lfl7ll . In the studies presented here, we rely 
only on charged particle tracking in the main drift cham- 
ber (MDC) and energy deposition in the electromagnetic 
calorimeter (EMC). 

The GEANT4-based Monte Carlo (MC) simulation soft- 
ware, which includes the geometric description of the BE- 
SIII detector and the detector response, is used to optimize 
the event selection criteria, determine the detection efficiency, 
and estimate backgrounds. For the signal process, we use a 
sample of e + e~ — > 7r + 7r~ J/tp MC events generated assum- 
ing the 7r + 7r~ J/tp is produced via F(4260) decays, and using 
the e + e~ — > tt + it~ J/tp cross sections measured by Belle H 
and BaBar The 7r + 7r~ J/ip substructure is modelled ac- 
cording to the experimentally observed Dalitz plot distribution 
presented in this analysis. ISR is simulated with KKMC ifljtl 
with a maximum energy of 435 MeV for the ISR photon, cor- 
responding to a 7r + 7r~ J/tp mass of 3.8 GeV/c 2 . 



For e + e~ — > ir + ir~ J/ip events, the J/ip candidate is re- 
constructed with lepton pairs (e + e~ or Since this de- 
cay results in a final state with four charged particles, we first 
select events with four good charged tracks with net charge 
zero. For each charged track, the polar angle in the MDC must 
satisfy | cos 9\ < 0.93, and the point of closest approach to the 
e + e _ interaction point must be within ±10 cm in the beam di- 
rection and within 1 cm in the plane perpendicular to the beam 
direction. Since pions and leptons are kinematically well sep- 
arated in this decay, charged tracks with momenta larger than 
1.0 GeV/c in the lab frame are assumed to be leptons, and 
the others are assumed to be pions. We use the energy de- 
posited in the EMC to separate electrons from muons. For 
muon candidates, the deposited energy in the EMC should be 
less than 0.35 GeV; while for electrons, it should be larger 
than 1 . 1 GeV. The efficiencies of these requirements are de- 
termined from MC simulation to be above 99% in the EMC 
sensitive region. 

In order to reject radiative Bhabha and radiative dimuon 
(7e + e~/7/i + /i~) backgrounds associated with a photon- 
conversion, the cosine of the opening angle of the pion can- 
didates, which are true e + e~ pairs in the case of background, 
is required to be less than 0.98. In the e + e~ mode, the same 
requirement is imposed on the 7r ± e T opening angles. This 
restriction removes less than 1% of the signal events. 

The lepton pair and the two pions are subjected to a 
four-constraint (4C) kinematic fit to the total initial four- 
momentum of the colliding beams in order to improve the mo- 
mentum resolution and reduce the background. The \ 2 of the 
kinematic fit is required to be less than 60. 

After imposing these selection criteria, the invariant mass 
distributions of the lepton pairs are shown in Fig. Q] A clear 
J/tp signal is observed in both the e + e~ and modes. 



There are still remaining e H 



7T + 7r 7r + 7r , and other 



QED backgrounds, but these can be estimated using the events 
in the J/tp mass sideband. The final selection efficiency is 
53.8 ± 0.3% for \i+ events and 38.4 ± 0.3% for e+e~ 
events, where the errors are from the statistics of the MC 
sample. The main factors affecting the detection efficiencies 
include the detector acceptances for four charged tracks and 
the requirement on the quality of the kinematic fit adopted. 
The lower efficiency for e + e~ events is due to final-state- 
radiation (FSR), bremsstrahlung energy loss of e + e~ pairs 
and the EMC deposit energy requirement. 

To extract the number of 7r + 7r~ J / tp signal events, invariant 
mass distributions of the lepton pairs are fit using the sum of 
two Gaussian functions with a linear background term. The 
fits yield M{J/tp) = 3098.4 ± 0.2 MeV/c 2 with 882 ± 33 
signal events in the mode; and M(J/tp) = 3097.9 ± 

0.3 MeV/c 2 with 595 ± 28 signal events in the e + e~ mode. 
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FIG. 1: The distributions of M(fj, (i ) (left panel) and M(e + e ) (right panel) after performing a 4C kinematic fit and imposing all selection 
criteria. Dots with error bars are data and the curves are the best fit described in the text. 



Here the errors are statistical only. The mass resolution is 
3.7 MeV/c 2 in the mode and 4.0 MeV/c 2 in the e + e~ 

mode. 

The Born cross section is determined from the relation 



N 



—0 , where N is the number of signal events 
from the fit; £ mt is the integrated luminosity; e is the selection 
efficiency obtained from a MC simulation; B is the branching 
fraction of J/tp —> £ + £~; and (1+6) is the radiative correction 
factor, which is 0.818 according to a QED calculation lfl9ll . 
The measured Born cross section for e + e~ — » ir + Tr~ J/tp is 
(64.4 ± 2.4) pb in the fi+ mode and (60.7 ± 2.9) pb in 
the e + e~ mode. The combined measurement is a B (e + e~ — > 
n+TT-J/tp) = (62.9 ± 1.9) pb. 

Systematic errors in the cross section measurement come 
from the luminosity measurement, tracking efficiency, kine- 
matic fit, background estimation, dilepton branching fractions 
of the J/tp, and F(4260) decay dynamics. 

The integrated luminosity of this data sample was measured 
using large angle Bhabha events, and has an estimated uncer- 
tainty of 1.0%. The tracking efficiency uncertainty is esti- 
mated to be 1% for each track from a study of the control 
samples J/tp — > ir + ir~ir° and -0(3686) — > tt + tt~ J/tp. Since 
the luminosity is measured using Bhabha events, the tracking 
efficiency uncertainty of high momentum lepton pairs partly 
cancels in the calculation of the 7r + 7r - J/tp cross section. To 
be conservative, we take 4% for both the e + e~ and /i + /x~ 
modes. 

The uncertainty from the kinematic fit comes from the 
inconsistency between the data and MC simulation of the 
track helix parameters. Following the procedure described in 
Ref. i20ll . we take the difference between the efficiencies with 
and without the helix parameter correction as the systematic 
error, which is 2.2% in the mode and 2.3% in the e + e~ 

mode. 

Uncertainties due to the choice of background shape and fit 
range are estimated by varying the background function from 
linear to a second-order polynomial and by extending the fit 



range. 

Uncertainties in the F(4260) resonance parameters and 
possible distortions of the Y (4260) line shape introduce small 
systematic uncertainties in the radiative correction factor and 
the efficiency. This is estimated using the different line shapes 
measured by Belle and BaBar 10]. The difference in 
(1 + S) ■ e is 0.6% in both the e + e~ and modes, and 

this is taken as a systematic error. 

We use the observed Dalitz plot to generate F(4260) — > 
7r + 7r~ J/tp events. To cover possible modelling inaccuracies, 
we conservatively take the difference between the efficiency 
using this model and the efficiency using a phase space model 
as a systematic error. The error is 3.1% in both the fi + [i~ and 
the e + e~ modes. 

The uncertainty in B(J/ip — > i + £~) is 1% 12 111 . The trigger 
simulation, the event start time determination, and the final- 
state-radiation simulation are well-understood; the total sys- 
tematic error due to these sources is estimated to be less than 
1%. 

Assuming all of the sources are independent, the total sys- 
tematic error in the ir + ir~ J/tp cross section measurement is 
determined to be 5.9% for the mode and 6.8% for the 

e + e~ mode. Taking the correlations in errors between the two 
modes into account, the combined systematic error is slightly 
less than 5.9%. 

Intermediate states are studied by examining the Dalitz plot 
of the selected tt + it~ J/tp candidate events. The J/tp signal 
is selected using 3.08 < M (£+£-) < 3.12 GeV/c 2 and the 
sideband using 3.00 < M(£ + l~) < 3.06 GeV/c 2 or 3.14 < 
M(£ + Z~) < 3.20 GeV/c 2 , which is three times the size of 
the signal region. In total, a sample of 1595 ir + ir~ J/tp events 
with a purity of 90% is obtained. 

Figure |2] shows the Dalitz plot of events in the J/tp signal 
region, where there are structures in the ir + 7r~ system and ev- 
idence for an exotic charmoniumlike structure in the n^J/ip 
system. The inset shows background events from J/tp mass 
sidebands (not normalized), where no obvious structures are 
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FIG. 2: Dalitz distributions of M 2 (tt + tt~) vs. M 2 (% + J/ip) for 
selected e + e~ — > 7r + 7r~ J/-0 events in the J/i/> signal region. The 
inset shows background events from the J/ip mass sidebands (not 
normalized). 



observed. 

Figure [3] shows the projections of the M(ir + J/ip), 
M \n~ J I 'ip) , and M(7r + 7r~) distributions for the signal 
events, as well as the background events estimated from nor- 
malized J/ip mass sidebands. In the 7r ± J/ip mass spectrum, 
there is a significant peak at around 3.9 GeV/c 2 (referred to 
as the Z c (3900) hereafter). The wider peak at low mass is a 
reflection of the Z c (3900) as indicated from MC simulation, 
and shown in Fig. [3] Similar structures are observed in the 
e + e~ and separated samples. 

The 7r + 7r~ mass spectrum shows non-trivial structure. To 
test the possible effects of dynamics in the n + ir~ mass spec- 
trum on the 7r ± J / ip projection, we develop a parameterization 
for the 7r + 7r~ mass spectrum that includes a /o(980), er(500), 
and a non-resonant amplitude. An MC sample generated with 
this parameterization adequately describes the 7r + 7r~ spec- 
trum, as shown in Fig. [3] but does not generate any peak- 
ing structure in the tt J/ip projection consistent with the 
Z c (3900). We have also tested _D-wave ir + ir~ amplitudes, 
which are not apparent in the data, and they, also, do not gen- 
erate peaks in the 7r ± J/tp spectrum. 

An unbinned maximum likelihood fit is applied to the 
distribution of M max (7r ± J/ip), the larger one of the two 
mass combinations M(ir + J/ip) and M (tt~ J/ip) in each 
event. The signal shape is parameterized as an S-wave Breit- 
Wigner (BW) function convolved with a Gaussian with a mass 
resolution fixed at the MC simulated value (4.2 MeV/c 2 ). The 
phase space factorp g is considered in the partial width, where 
p is the Z c (3900) momentum in the Y (4260) CM frame and q 
is the J/ip momentum in the Z c (3900) CM frame. The back- 
ground shape is parameterized as a/(x — 3.6) b + c+dx, where 
a, 6, c, and d are free parameters and x = M max (ir J /ip)- 
The efficiency curve is considered in the fit and the possi- 
ble interference between the signal and background is ne- 
glected. Figure|4]shows the fit results; the fit yields a mass of 
(3899.0 ± 3.6) MeV/c 2 , and a width of (46 ± 10) MeV. The 



goodness-of-the-fit is found to be \ 2 /ndj = 32.6/37 = 0.9. 

The number of Z c (3900) events is determined to be 
A r (Z c (3900) ± ) = 307±48. The production ratio is calculated 
to be R = ^^0^^fIM = (21.5±3.3)%, 
where the efficiency correction has been applied. The statisti- 
cal significance is calculated by comparing the fit likelihoods 
with and without the signal. Besides the nominal fit, the fit is 
also performed by changing the fit range, the signal shape, or 
the background shape. In all cases, the significance is found 
to be greater than 8a. 

Fitting the M(ir + J/ip) and M(n~J/ip) distributions sep- 
arately, one obtains masses, widths, and production rates of 
the Z c (3900)+ and Z c (3900)~ that agree with each other 
within statistical errors. Dividing the sample into two dif- 
ferent A/(7r + 7r~) regions (below and above AI 2 (tt + it~) = 
0.7 GeV 2 /c 4 ) allows us to check the robustness of the 
Z c (3900) signal in the presence of two different sets of inter- 
fering 7r + 7r~ J /ip amplitudes. In both samples, the Z c (3900) 
is significant and the observed mass can shift by as much 
as 14 ± 5 MeV/c 2 from the nominal fit, and the width can 
shift by 20 ± 11 McV. We attribute the systematic shifts 
in mass and width to interference between the Z c (3900)7r 
and (TT + Tr~)J/ip amplitudes. In fitting the n^J/ip projec- 
tion of the Dalitz plot, our analysis averages over the en- 
tire 7r + 7r~ spectrum, and our measurement of the Z c (3900) 
mass, width, and production fraction neglects interference 
with other tt + it~ J/ip amplitudes. 

The systematic errors for the resonance parameters of the 
Z c (3900) come from the mass calibration, parametrization 
of the signal and background shapes, and the mass resolu- 
tion. The uncertainty from the mass calibration can be esti- 
mated using the difference between the measured and known 
J/ip masses (reconstructed from e + e~ and and D° 

masses (reconstructed from K~tt + ). The differences are 
(1.4 ± 0.2) MeV/c 2 and -(0.7 ± 0.2) MeV/c 2 , respectively. 
Since our signal topology has one low momentum pion, as 
in D° decay, and a pair of high momentum tracks from the 
J / ip decay, we assume these differences added in quadrature 
is the systematic error of the Z c (3900) mass measurement 
due to tracking. Doing a fit by assuming a P-wave between 
the Z c (3900) and the ir, and between the J/ip and tt in the 
Z c (3900) system, yields a mass difference of 2.1 MeV/c 2 , 
a width difference of 3.7 MeV, and production ratio differ- 
ence of 2.6% absolute. Assuming the Z c (3900) couples 
strongly with DD* results in an energy dependence of the 
total width J22j], and the fit yields a difference of 2.1 MeV/c 2 
for mass, 15.4 MeV for width, and no change for the pro- 
duction ratio. We estimate the uncertainty due to the back- 
ground shape by changing to a third-order polynomial or a 
phase space shape, varying the fit range, and varying the re- 
quirements on the x 2 of the kinematic fit. We find differences 
of 3.5 MeV/c 2 for mass, 12. 1 MeV for width, and 7. 1 % abso- 
lute for the production ratio. Uncertainties due to the mass 
resolution are estimated by increasing the resolution deter- 
mined by MC simulations by 16%, which is the difference 
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FIG. 3: One dimensional projections of the M(tv + J/tp), M(ir~ J/tp), and A'f(7r + 7r~) invariant mass distributions in e + e~ — ► tt + -k~ J/tp 
for data in the J/tp signal region (dots with error bars), data in the J/tp sideband region (shaded histograms), and MC simulation results from 
<j(500), /o(980) and non-resonant tt + -k~ amplitudes (red dot-dashed histograms). The pink blank histograms show a MC simulation of the 
Z c (3900) signal with arbitrary normalization. 
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quark. Similar studies were performed in B decays, with un- 
confirmed structures reported in the 7r ± i/;(3686) and it Xc\ 
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FIG. 4: Fit to the M ma x(7T J/tp) distribution as described in the 
text. Dots with error bars are data; the red solid curve shows the total 
fit, and the blue dotted curve the background from the fit; the red dot- 
dashed histogram shows the result of a phase space MC simulation; 
and the green shaded histogram shows the normalized J/tp sideband 
events. 



between the MC and measured mass resolutions of the J / ip 
and D° signals. We find the difference is 1.0 MeV in the 
width, and 0.2% absolute in the production ratio, which are 
taken as the systematic errors. Assuming all the sources of 
systematic uncertainty are independent, the total systematic 
error is 4.9 MeV/c 2 for mass, 20 MeV for width and 7.5% for 
the production ratio. 

In Summary, we have studied e + e~ — > it + tt~ J/tp at a 
CM energy of 4.26 GeV. The cross section is measured to 
be (62.9 ± 1.9 ± 3.7) pb, which agrees with the existing re- 
sults from the BaBar Jl, Belle Jj, and CLEO H experi- 
ments. In addition, a structure with a mass of (3899.0 ± 3.6 ± 
4.9) MeV/c 2 and a width of (46 ± 10 ± 20) MeV is observed 
in the w J/tp mass spectrum. This structure couples to char- 
monium and has an electric charge, which is suggestive of a 
state containing more quarks than just a charm and anti-charm 



systems 123142611 . It is also noted that model-dependent calcu- 
lations exist that attempt to explain the charged bottomonium- 
like structures which may also apply to the charmoniumlike 
structures, and there were model predictions of charmonium- 
like structures near the DD* and D*D* thresholds 027I1. 
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